Singular conductance of a spin 1 quantum dot 
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We interpret the recent observation of a zero-bias anomaly in spin-f quantum dots in terms of an 
underscreened Kondo effect. Although a spin-f quantum dots are expected to undergo a two-stage 
quenching effect, in practice the log normal distribution of Kondo temperatures leads to a broad 
temperature region dominated by underscreened Kondo physics. General arguments, based on the 
asymptotic decoupling between the partially screened moment and the leads, predict a singular 
temperature and voltage dependence of the conductance G and differential conductance g, resulting 
in dg/dT ~ 1/T and dG/dV ~ 1/V. Using a Schwinger boson approach, we show how these 
qualitative expectations are borne out in a detailed many body calculation. 
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Single-electron transistors (SETs) offer the intriguing 
opportunity to probe and explore classes of strongly cor- 
related electron behavior associated with the Kondo ef- 
fect that are difficult to access in bulk materials 0,0 00 
5] . The possibility of observing a break-down in Landau 
Fermi liquid behavior that accompanies the overscreened 
two-channel Kondo effect in quantum dots has been a 
subject of particular recent interest^, 0]. In this paper, 
we propose that singular deviations from Landau Fermi 
liquid behavior associated with the underscreened Kondo 
effect, hitherto unobserved in bulk materials, will develop 
in conventional quantum dots with even numbers of elec- 
trons and a triplet ground-state 0, 0. These devi- 
ations from conventional Fermi liquid behavior are pre- 
dicted to lead to singular voltage, field and temperature 
dependences in the conductance. 

The Kondo effect in quantum dots with odd numbers 
of electrons, predicted more than fifteen years ago, 00 
is now well-established by experiment HI ■ Subsequent 
observations have shown that zero-bias anomalies associ- 
ated with a Kondo effect can also occur in quantum dots 
with even occupancies, where Hund's coupling between 
the electrons can lead to novel degeneracies, through the 
formation of higher spin states, or the accidental de- 
generacy of singlet and triplet states. Zero-bias anoma- 
lies in integer spin quantum dots were first reported by 
Schmid et al.Q. Sasaki et al@ later discovered a zero- 
bias anomaly in even electron quantum dots, associated 
with the degeneracy point between singlet and triplet 
states, tuned by a small magnetic field. Most recently, 
Kogan et al|I3* have shown that the singlet-triplet exci- 
tation energy in lateral quantum dots can be tuned by 
the gate voltage, explicitly demonstrating that the zero 
bias anomaly develops once the triplet state drops below 
the singlet configuration. 

Pustilnik and Glazman|f l| have analyzed the low- 
temperature Fermi liquid physics of higher spin quantum 
dots. Their analysis shows that lateral quantum dot in 
a triplet configuration develops two screening channels 



which fully screen the local moment at the lowest tem- 
peratures. Using the Landauer formula, they deduce the 
conductance G of the Fermi liquid which develops to be 

G=^-sin 2 (£ 1 -£ 2 ) (f) 

where 5\ and 62 are the scattering phase shifts of the two 
screening channels. According to this line of reasoning, 
the development of a unitary phase shift in each channel, 
Si =82 = 7r/2 leads to a complete suppression of the zero 
bias anomaly in a triplet quantum dot JL2| . Why then, are 
zero-bias anomalies, with near unitary conductance seen 
in triplet quantum dots? 

fn this paper we propose an interpretation of this un- 
expected behavior in terms of an underscreened Kondo 
effect. Our key observation is that the antiferromagnetic 
Kondo coupling constants J\ (A = 1,2) associated with 
the two screening channels in a triplet quantum dot will 
generally be distributed independently. Since the Kondo 

temperature depends exponentially on the coupling con- 
_ 1 

stant Tk\ = D\fJ\pe ■'>■'' , a normal distribution of the 
coupling constants will drive a log -normal distribution in 
the two Kondo temperatures [fjj, with the potential to 
generate exponentially large separations in the relative 
magnitude of the Kondo temperatures of each channel. 
If we assume that ln(Txi/T^ 2 ) = jjr — 77^ >> 1, then 
over the exponentially broad temperature range given by 
log(Tftri) >> logT >> log(7if 2 ), the underlying physics 
is that of a one channel spin-f Kondo model, in which 
the spin is partially screened to a spin 1/2. 

From this perspective, triplet dots with a large zero 
bias anomaly are those where the Kondo coupling con- 
stants of the two channels are severely mis-matched, giv- 
ing rise to decades of behavior dominated by the under- 
screened Kondo effect in a single channel. Previous work, 
both analytic |ll| and numerical has focussed on 

the equilibrium behavior of triplet quantum dots with 
Kondo temperatures of comparable magnitude. We now 
examine the singular consequences of a wide separation 
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between these two scales in both finite temperature and 
finite voltage properties. 

In the underscreened spin-1 Kondo effect, the residual 
spin-1/2 moment is ferromagnetically coupled to leads, 
with a coupling that scales logarithmically slowly to 
zero|l5j. The ground-state which develops is a "singular 
Fermi liquid" , in which the electrons do behave as Lan- 
dau quasiparticles which are elastically scattered with 
unitary phase shift, but where, on the other hand, a loga- 
rithmically decaying coupling generates a singular energy 
dependence in the scattering phase shift and a di verg ence 
in the resulting quasiparticle density of states flo. fl7t fl8|| . 
The Hamiltonian for the underscreened quantum dot is 



H = ^ e k cl Xa c kX a 



(2) 



where ipa = J2k ac kLa + @CkRa denotes the linear combi- 
nation of right and left channels that couples to the dom- 
inant screening channel. Much is known about the equi- 
librium physics of this model. At low temperatures, the 
spin is partially screened from spin S to spin S — (1/2). 
The residual moment is ferromagnetically coupled to the 
conduction sea, with a residual coupling that slowly flows 
to weak coupling according to 



Jp(A) 
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(3) 



where A <~ max(T, (AbB) is the characteristic cut-off en- 
ergy scale, provided in equilibrium, by the temperature 
or magnetic field. At low energies and temperatures, the 
partially screened magnetic moment scatters electrons 
elastically, with a unitary phase shift, however the cou- 
pling to the residual spin (S — ^) gives rise to a singu- 
lar energy dependence of the scattering phase shift. The 
low energy scattering phase shift can be directly deduced 
from the Bethe Ansatz, and has the asymptotic form 
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The logarithmic term on the right hand side is pro- 
duced by the residual coupling between the electrons 
and the partially screened moment. While the elec- 
trons at the Fermi energy scatter elastically off the lo- 
cal moment with unitary scattering phase shift, as in 
a Fermi liquid, the logarithmically singular dependence 
of the phase shift leads to a divergent density of states, 



7T du 



t^t, which means that we can not as- 
sociate this state with a bona-fide Landau Fermi liquid. 
For this reason, the ground-state of the underscreened 
Kondo model has recently been called a "singular Fermi 
liquid" [13. 



These singular features of the underscreened Kondo ef- 
fect are expected to manifest themselves in the properties 
of the triplet quantum dot. For example, we expect the 
low-field conductance to follow the simple relation 



G(B) 



2e2 • 2 X / d n 2 e 2 



1 



16 ln 2 ( 



(5) 



for B « Tk- This relationship was previously obtained 
by other means from the Tki — > of the two-channel 
model|ll|. Notice that the field derivative of the conduc- 
tance diverges as ^ oc 1/ (bIii 3 (Tk/B)^ at low fields. 

The prediction of the finite temperature, and finite volt- 
age conductance can not be made exactly, however we 
expect the above form to hold, for the differential con- 
ductance at finite temperature or voltage, with an appro- 
priate replacement of cut-offs, namely 
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(6) 



and dG/dT ~ l/max(T, V). 

To model this behavior in more detail it is useful to 
consider a simplified model of the quantum dot in which 
the Hund's coupling is taken to be infinite. In this limit, 
the states of the quantum dot can be described using a 
Schwinger boson representation as 



\d\a) 



\d 2 ,M)=b^ M _ a \0), (7) 



Written in this representation the model becomes 

H = H +tY / fylxK + b\ X ^a\ + E dX h (8) 

a 

subject to the constraint + x'x = 2. 

To develop a controlled many body treatment of this 
Hamiltonian, we use a large-A^ expansion, extending the 
number of spin components a from two to N. To pre- 
serve a finite scattering phase shift as N — > oo, we intro- 
duce K = kN bosonic "replicas" , where k is fixed. With 
this device we obtain a (dynamical) mean field theory 
with scattering phase shift S = irk and the qualitatively 
correct logarithmic energy dependences 01 • The Hamil- 
tonian used in the large N expansion is then 



H = H Q 



kN 



xk/*-(9) 



In the large N limit, there are two self-consistent non- 
crossing approximations to the Dyson equations for the 
self-energies of the conduction electrons and x fermions 
(Fig. 1). Here we sketch the main elements of the deriva- 
tion. As in the corresponding equilibrium calculation^^], 
the boson behaves as a sharp excitation in the large N 
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limit, with an average occupancy {n^^) — nb/N. From 
the Dyson equations we obtain sets of self-consistent in- 
tegral equations for both the retarded and Keldysh self- 



energies. The explicit expressions for the the retarded 
self-energies S^. and H C R of the slave fermion (x) and the 
conduction electrons (c) are 



t UfeG^CA — uj) — t 
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(10) 



Here Jr(lu) = [u - E d - £%H] _1 and G%(u) = 
[(— irp)^ 1 — are the retarded propagators for 

the x fermions and conduction electrons. 
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FIG. 1: The non-crossing approximation for the self-energies 
of conduction electrons and \ fermions. The solid line denotes 
the Larkin Ovchinnikov matrix propagator for the conduction 
electrons. The dashed line denotes the corresponding Green's 
function of the auxiliary (%) fermions and the wavy line is the 
bosonic propagator. Thin lines denote the bare propagator 
and full lines the dressed propagator. Each vertex corresponds 
to the factor . 

The ratio of the Keldysh to the retarded self-energies 
self-consistently determines the fermion distribution 
functions. We can summarize the results of our calcu- 
lation of the Keldysh self-energies by providing the dis- 
tribution functions that they generate. The distribution 
function of the conduction electrons is the average 



/c = ^[/iH + /*(*)] 



(ii) 



where !l,r{oj) = l/( e l3( -^ eV / 2 + 1) is the equilibrium 
distribution function in the left/right-hand lead. The 
distribution function of the auxiliary fermion is 



n b 



(12) 



where = l/(e^ A + 1) determines A. This relationship 
can be simply understood as the result of detailed balance 
between rate of the decay processes c — > b + \ an d b + 




FIG. 2: Imaginary part of dot T-matrix for a variety of 
voltages for the case k = 0.4. As the voltage is increased, the 
singular central peak splits into two components. 



X —> c, and it reverts to the equilibrium Fermi Dirac 
distribution in the limit V — > 0. 

From these results, we compute the temperature and 
voltage dependent current, given by 0] 



2 V 



I(V,T) = N—p / dw 
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where £r(w) = T, r (uj)/[1 — itt pT, R (u))] is the scattering 
t-matrix. 

We have solved these equations numerically, and the 
key results are shown in Figs 2-4. Fig. 2. shows the volt- 
age dependent t-matrix at zero temperature. At zero 
voltage, the t-matrix contains a logarithmic singular- 
ity noted in previous work[ll [3 

which splits into two 
peaks at a finite voltage. In our calculation, the split 
Kondo resonance retains its singular structure, although 
this is most likely an artifact of taking a limit where 
the bosons behave as a sharp excitation. In Fig. 3., 
we show the temperature dependent conductance. The 
temperature-dependent deviations from unitary conduc- 
tance are determined by the logarithmic singularity in 
the phase shift, and in our calculation, these are propor- 
tional to l/ln(Tx/max(eV, T)). In the Schwinger boson 
approach, the number of bound bosons in the Kondo sin- 
glet never exceeds N/2 and the region K > N/2 does not 
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onance is predicted to give rise to a 1 /T divergence in the 
temperature dependence of the differential conductance, 
and a 1/V divergence in the second derivative of the 
voltage dependent current d 2 I/dV 2 . These ideas have 
been developed qualitatively and illustrated within an 
integral equation treatment of the underscreened Kondo 
model. Experimental observation of these singular fea- 
tures would constitute a first realization of the under- 
screened Kondo effect. 
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00ER45790 (PC). 



FIG. 3: Temperature dependence of the differential conduc- 
tance, normalized with respect to gu = N^-sin 2 (irk) for the 
representative case k = 0.4. Insert shows the 1/T divergence 
of the derivative dg/dT. 

describe an underscreened Kondo model. Consequently, 
we are limited to static phase shifts 6 — tt(K/N) < ir/2, 
so the strictly particle- hole symmetric case S = it/2 is 
outside the limits of our approach. Nevertheless, our 
numerical results do capture the expected singularities. 
Fig. 3. shows the singular form of the temperature de- 
pendent conductance, with singular 1/T divergence in 
dg/dT. Finally, Fig. 4. shows the voltage dependence of 
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FIG. 4: Voltage dependent conductance G(V) = I(V)/V for 
the case k = 0.4. Insert: dG(V) / dV showing 1/V divergence. 

the conductance, which has a similar logarithmic singu- 
larity at low voltage. 

In summary, we have proposed that the monotonically 
increasing conductance observed as the temperature is 
lowered in triplet quantum dots is associated with an 
underscreened Kondo effect. The singular energy and 
temperature dependence associated with the Kondo res- 
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